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Abstract 

Thermal history of the Universe between inflation and big-bang nucleosynthesis has not yet 
been revealed observationally. It will be probed by the detection of primordial gravitational waves 
generated during inflation, which contain information on the reheating temperature as well as 
the equation of state of the Universe after inflation. Based on Fisher information formalism, we 
examine how accurately the tensor-to-scalar ratio and reheating temperature after inflation can be 
simultaneously determined with space-based gravitational wave detectors such as the DECI-hertz 
Interferometer Gravitational- wave Observatory (DECIGO) and the Big-Bang Observer (BBO). We 
show that the reheating temperature is best determined if it is around lO^GeV for tensor-to-scalar 
ratio of around 0.1, and explore the detectable parameter space. We also find that equation of state 
of the early Universe can be also determined accurately enough to distinguish different equation- 
of-state parameters if the infiationary gravitational waves are successfully detected. Thus future 
gravitational wave detectors provide a unique and promising opportunity to reveal the thermal 
history of the Universe around lO^GeV. 
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I. INTRODUCTION 



The standard paradigm in cosmology is well established and explains cosmological evolu- 
tion after the big-bang nucleosynthesis (BBN). An epoch of primordial inflation at the very 
early Universe now also becomes an essential part of the standard paradigm, and is supported 
by accurate measurements of the cosmic microwave background (CMB) anisotropics. The 
nearly scale-invariant power spectrum of the curvature perturbation, as well as its Gaussian 
nature, strongly indicates the inflation as a source of the density fluctuation of the Universe 

Q. 

However, thermal history between inflation and BBN has not yet revealed observationally. 
There must be an epoch of reheating after inflation, where the inflaton decays and the 
radiation-dominated universe starts. The process of reheating crucially depends on the 
properties of the inflaton: its mass, potential and couplings to standard model particles. 
Before the inflaton decays, inflaton oscillates around its potential minimum. In the simplest 
class of inflation models, the inflaton oscillation behaves as matter after inflation ends. 
But the detail of the inflaton oscillation, or the equation of state of the Universe during 
the inflaton oscillation depends on the nature of the inflaton. Therefore, it is extremely 
important to probe the thermal history after inflation. It would determine the nature of the 
inflaton and might be a direct hint for the underlying high-energy theory such as supergravity 
or string theory. 

Then, how can we probe the very early Universe observationally? Perhaps gravitational 
waves may be a unique signal that directly carries information on the thermal history of 
these early epochs, since the Universe is transparent to gravitational waves up to the Planck 
epoch in principle. Fortunately, there is a promising source of the gravitational wave: pri- 
mordial gravitational waves generated by inflation [2,]. Quantum fluctuations during inflation 
are expanded and frozen outside the Hubble horizon. Tensor fluctuations generated in this 
way can be observed as gravitational waves. These gravitational waves form a stochastic 
background having a very wide range of wavelength from the present horizon scale to the 
terrestrial one. If the inflationary energy scale is sufficiently large, inflationary gravitational 
waves will be detected through on-going or future CMB B-mode polarization measurements 
0,0] and/or space laser interferometers js, 6|. Proposed future gravitational wave measure- 
ments such as the DECI-hertz Interferometer Gravitational-wave Observatory (DECIGO) 
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B, B and the Big-Bang Observer (BBO) fl Q aim at the detection of the inflationary 
gravitational wave background around 0.1-1 Hz as well as the determination of the nature 
of the present accelerated expansion of the Universe with binary sources as standard sirens 
lo|-|l2|. 

Interestingly, the spectrum of the inflationa ry g ravitational wave background directly 



reflects the expansion history of the Universe 13l4l6|. This is because the gravitational wave 



amplitude frozen at super-horizon scale begins to damp inversely proportional to the scale 
factor of the Universe when the corresponding mode enters the horizon. As we have already 
discussed, the matter-dominated reheating period is followed by the radiation-dominated 
era. Therefore, the spectrum shows a knee-like feature at the frequency corresponding to the 
comoving horizon scale at the end of reheating, which is given by /r ^ 0.26 Hz(Tr/10^ GeV) 
(see Eq. f lT3|) below) where Tr is the temperature of the Universe at the end of reheating 
defined by 

Tr = , yrvi^, (1) 

with being the inflaton decay rate and Mp being the reduced Planck scale. If future grav- 
itational wave detectors see the knee shape at this frequency, it will be a direct measurement 
of the reheating temperature. 

gravitational wave background with 



17Hl9l|. The signal-to- noise ratio is 



The prospects for direct detection of the inflationary 
the effect of reheating have been investigated in Refs. 
calculated taking into account the effect that reheating induces a change of the frequency 
dependence in the spectrum, and the detectability of the gravitational wave background 
is discussed in the previous works. However, direct detection of the gravitational wave 
amplitude does not necessarily mean accurate determination of the reheating temperature. 
In order to estimate the detectability of the reheating signature more realistically, we need 
to examine if the experiment can distinguish the change of the characteristic frequency 
dependence induced by reheating. For that purpose, we employ the Fisher matrix formalism, 
which enables us to evaluate the sensitivity to the shape of the spectrum with properly taking 
into account parameter degeneracies. 

As we have mentioned, the change of the frequency dependence is attributed to the 
fact that reheating causes the change of the Hubble expansion rate. Thus, the inflationary 
gravitational wave background is useful not only determining the reheating temperature but 
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also constraining the expansion history of the very early Universe. We evaluate the ability of 
future direct detection experiments to determine the equation of state of the early Universe, 



keeping in mind that many applications on cosmological mode 
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s are expected with direct 



23|. 



detection of the inflationary gravitational wave background 

In this paper, we study prospects for determination of the reheating temperature, as well 
as the equation of state of the early Universe, at future gravitational wave detectors based 
on the Fisher matrix method. In Sec. Ull basic formulations are provided. In Sec. Illlt we 
derive the accuracy of the reheating temperature measurement with future gravitational 
wave detectors. In Sec. HV] we study how accurately the equation of state of the early 
Universe will be determined in the same setup. The final section is devoted to conclusions. 
The fiducial cosmological parameters are taken to be the maximum likelihood values from 
the combined analysis of the WMAP 7-year, baryon acoustic oscillation, and supernova 
data with a fiat ACDM universe jl]: matter density fim^^ = 0.1344, amplitude of curvature 
perturbation A^ pj,jj^ = 2.45 x 10~^ and the Hubble parameter h = 0.702. Throughout the 
paper we choose the natural units, c = h = 1. 



II. FISHER MATRIX METHOD 



A. Definitions of gravitational wave observables 

Before discussing the Fisher matrix approach, let us first summarize the observable 
quantities of the gravitational wave detectors. Gravitational waves in the expanding Uni- 
verse are described as tensor perturbations in the Friedmann-Robertson- Walker metric, 
ds^ = a'^{t)[—dT'^ + {6ij + hij)dx'^dx^], with a{t) being the scale factor of the Universe. The 
conformal time r is defined as dr = dt/a{t), and the subscript i and j denote spatial indices 
which run over 1, 2, 3 or x, y, z. The tensor perturbation hij satisfies the transverse-traceless 
conditions, d^hij = h] = 0. It is expanded into its Fourier series as 

^^^■(^'-)= E / 7^4■(k)/^^We^'^■^ (2) 

where the polarization tensors ep^ satisfy symmetric and transverse-traceless conditions 
and are normalized as j ^iji^ij)* = 26^'^'. Then the observable intensity of a stochastic 
gravitational wave background is characterized by its density parameter per logarithmic 
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wavenumber (or frequency) as 

1 rfPGW ^ / \ 2 Z.3 



pcdlnk 12 \aH J tt^ 

where pcw denotes the energy density of the gravitational waves, pcw = {{drhijY + 
(V/ijj)^)/(647rG'a^), = ^H'^/SnG is the critical density of the Universe, and the Hub- 
ble parameter is defined a.s H = {da/dt)/a. One can also express such a statistical quantity 
in terms the power spectrum, 

A 

B. Fisher matrix for gravitational wave measurement 

Here, we briefly review and discuss the Fisher information matrix approach, focusing 
on the future gravitational wave detectors. The Fisher matrix is a powerful method to 
theoretically forecast the constraining power on parameters of interest for a given survey. 



and is commonly used in observational cosmology (see e.g., [2J-|26| ). An essential but unique 
assumption in the formalism is a Gaussian likelihood, and the Fisher matrix is defined by 
the second derivative of the likelihood around its maximum (or fiducial parameters) with 
respect to the parameters of interest. Then the inverse of the Fisher matrix provides a lower 
bound on the covariance matrix via the Cramer-Rao bound, and is regarded as the best 
achievable accuracy for the parameters. 

The Fisher matrix generally depends on the covariance matrix of signals, i.e., the noise 
properties of the survey configurations. In the case of the stochastic gravitational wave back- 
ground, direct detection can be attempted by cross- correlating the output signals between 
two detectors For the cross-correlation analysis, the Fisher information matrix is given 



by |28| 



/3^y >^ ff- Jlij{f)\%MGw{m^nay,if) 

with / = k/{2TT) and Hq being the frequency of the gravitational wave and the present 
Hubble parameter. We choose a lower cutoff of /cut = O.lHz, below which the signal may 



be contaminated by noise from cosmological white dwarf binaries 29|. As for /max, we set 
/max = oo unless otherwise stated, though the high-frequency range is limited by the noise 
spectrum which we will see in detail later. For a given survey, we need to assign three 



observational ingredients; a total observational Tobs, overlap reduction functions 'yij{f), and 
noise spectra Sr(.f)- These functional forms rely on the type of interferometry that we choose 



(for a review, see |30|). 

In the case of the Time-Delay Interferometry (TDl) that is expected to be adopted in 
BBO, the subscript / or J denotes the TDl channel output index {I, J = A,E,T). We 
compute the overlap reduction function '~fij{f) for the TDl data combinations with the 



method in Ref. j3l|- The noise transfer functions for the TDl variables are assumed as 32 1 
SAif) = SEif) = 8sin2(/72)[(2 + cos/^shot + 2(3 + 2cos/ + cos(2/))5accci], (5) 

5t(/) =2[l + 2cos/T[^shot + 4sin2(/72)5accei], (6) 

where Sghot and S'accei are the photon shot-noise and the proof-mass acceleration noise for 
the laser interferometers. The noise spectrum is designed so as to be rescaled with respect to 
the pivot frequency / = 2TcLf with L being the arm length. Unlike BBO, DECIGO would 
install a Fabry-Perot type interferometer, and hence its noise functions can be different from 
those of BBO discussed above. However, both DECIGO and BBO are basically designed to 
aim at the detection of the inflationary gravitational waves with similar frequency ranges 
around 0.1 — IHz. In addition, the purpose of this paper is to quantify the potential of future 
gravitational wave experiments for determination of model parameters. Based on these facts, 
we simply assume the TDl- type noise functions, Eqs. and ([HD, throughout the paper. 
Even with such a simple assumption, we do not believe that detailed configurations can 
significantly affect our conclusion. 

In table [U we summarize the survey parameters adopted in the Fisher matrix calculation. 
We consider two types of detectors. The first one is corresponding to the proposed DECIGO 
or BBO, quoted by "BBO/FP-DECIGO" , whose parameter values are taken from Ref. 33|. 
We also investigate the case of an ideal experiment whose sensitivity is limited only by 
quantum noises. We quote it as "Ultimate-DECIGO" and its parameters are taken from 
Ref. jsi]. 

In summary, the Fisher matrix can be calculated with Eq. ([5]), once the detector param- 
eters are assigned and theoretical predictions for ^gw is provided. Then the marginahzed 
la error is easily computed with the inverse of the Fisher matrix, 

^{P^) = ^/{^^^. (7) 
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*^shot 


*--'accel 


L 


Detectors 


[(L/km)-2Hz-i] 


[(27r//Hz)-4(2L/km)-2Hz-i] 


[km] 


BBO/FP-DECIGO 


2 X 10-^° 


9 X 10-40 


5 X 104 


Ultimate DECIGO 


9 X 10-4^ 


9 X 10-44 


5 X 104 



TABLE I: detector parameters 

The Fisher matrix is a product of the signal-to-noise-ratios, ficw/'S', and the derivative, 
9p. In ficW) and hence depends on the parameter response, namely the parameter de gene racy 



as well as the signal detectability. This makes a difference from the previous work 



III. DETERMINATION OF THE REHEATING TEMPERATURE 



18| 



As discussed in 
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19j, direct detection of the inflationary gravitational wave background 
has a potential to constrain or even to determine the reheating temperature via the char- 
acteristic frequency dependence. Inflationary gravitational wave spectrum has a frequency 
dependence of /"^ for modes which enter the horizon during a matter-dominated universe, 
and /° for modes which enter the horizon during a radiation-dominated universe. If the 
Universe behaves like a matter-dominated universe during reheating, the transition from 
reheating to the radiation domination is seen as a change of the frequency dependence of 
the spectrum. If this signature exists in the frequency band of direct detection sensitivity, 
0.1 — IHz, we may be able to determine the reheating temperature by measuring the knee- 
like feature where the frequency dependence changes from /"^ to /°. In this section, we 
apply the Fisher matrix analysis to investigate to what extent the future gravitational wave 
experiments can determine the reheating temperature. 

A. The gravitational wave spectrum in the presence of reheating 

The spectrum of the inflationary gravitational wave background is often expressed in 
terms of the initial tensor power spectrum ^\^^-^^JJt) and the transfer function Th{k), 



^G^ = Y2\-^) KpriAkmk). (8) 

In a single field slow-roll inflation, the tensor-to-scalar ratio r = A^ pj,;jj^(fco)/A^ pj,jj^(A;o) can 
be related to the tilt of the tensor mode spectrum ut = din Al^^^^^^kQ) / din k as r = — Sn^. 
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From this relation, the initial power spectrum can be written as 



^h,pnm{k) ^ rA^^ p,i^(A;o) exp 



r , k 
--In — + 

8 ko 



(9) 



where the pivot scale is taken as ko = 0.002 Mpc ^. ^ The effects of the cosmological 



evolution after inflation are all included in the transfer function, which is given as 18|, 121 1 
T'M = n?„ (*<5=)) f^)"' (^^f^X T^MTiM. (10) 

V 9*0 J \9*s{-Lm)J \ kTo J 

The subscript "0" denotes the present time and "in" denotes the time when the mode 
k crosses the horizon. The effective number of degrees of freedom at the end of reheat- 
ing is taken to be the sum of the standard model particles, 9*{Tb,) = (7*5 (Tr) = 106.75. 
The values at present are 9^0 = 3.36 and 9^so = 3.90. Here, tq is the present confor- 
mal time calculated assuming the Universe is matter dominated: tq = 2Hq^. The ef- 
fect of the cosmological constant is accounted for by the factor of Qm = 1 — ^A- In the 
limit of kTQ ^ 1, the spherical Bessel function ji{x) = (sinx — xcosa;)/x^ is replaced as 



jiikro) — i- cos(A;ro)/(A;ro) = l/{\/2kTQ). The first transfer function Ti(xcq) describes the 
change of the frequency dependence of the spectrum which arises from the change of the 
expansion rate of the Universe at the matter-radiation equality t = teq [sl , 



T^{xeq) = (1 + 1.57xeq + 3.42x,\), (11) 
where Xeq = k/k^q and k^q = = 7.1 x 10~^f2m^^ Mpc~^. The second transfer function 



T2{x^) corresponds to the change of the expansion rate at the end of reheating t = |18| . 



T^ixB.) = (1 - 0.32XR + 0.99xi)-\ (12) 

where xr = k/k^ and ^r ~ 1.7 x 10i^Mpc~^(^,,(TR)/106.75)i/''(TR/10^GeV). This can be 
rewritten in terms of frequency as 

^, n o«TT„ f 9*s{Tr)Y^^ f Tr 



^ The contribution from the higher order of ln(fc/fco) is sometimes non-neghgible, depending on the inflation 
model, and may cause a wrong estimate of the amplitude of the spectrum [35| . It can affect determination 
of the value of r, which is an important parameter to determine the gravitational wave amplitude, but is not 
crucial for determination of the reheating temperature, because the reheating temperature is determined 
basically by the characteristic frequency dependence of the spectrum. 
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FIG. 1: The spectra of the inflationary gravitational wave background for different values of 
the reheating temperature (thick black solid curves with Tr = 10^''''^ GeV from left to right). 
The tensor-to-scalar ratio is taken to be r = 0.1. For reference, the noise spectra for BBO/FP- 
DECIGO with 10-year observation (red solid) and for Ultimate-DECIGO with 3-year observation 
(blue dotted) are shown. The gray shaded region is not used in the Fisher analysis, since noises 
from white dwarf binaries may significantly contribute as systematic errors. 

which is the frequency where the change of the frequency dependence due to reheating arises. 
We show the spectra for different values of the reheating temperature in Fig. [H As clearly 
seen from the figure, the knee shape around fji can be observed with BBO/FP-DECIGO if 
the observational time is sufficiently long. Since there are no observable to probe reheating 
so far, the future gravitational wave experiments may provide us an unique opportunity to 
reveal the reheating of the Universe. 

B. Result 

Based on the theoretical prediction presented above, let us estimate the detectability 
of the reheating temperature using the Fisher matrix, which is calculated by substituting 
Eqs. (l8l)-(lT2l) into Eq. (jS]). We take r and Tr as free parameters, which correspond to the 
amplitude of the spectrum and the frequency of the reheating signature. 

In Fig. 121 we present an example of the expected future constraints, in which the fiducial 
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parameters are chosen as r = 0.1 and Tr = lO^GeV. Each elhpse represents the 2a error 
contours expected from 1, 3 and 10 years of observation with BBO/DECIGO. The error 
elhpse shrinks more for longer observations due to the fact that the signal-to-noise ratio 



scales as y^^jbs- As naturally expected, there is a degeneracy between Tr and r. This is 
simply because, as long as the frequency dependence of the spectrum is measured with a 
good accuracy we cannot distinguish the spectrum with larger Tr from that with smaller 
amplitude, i.e., with smaller r. 

An interesting and nontrivial question is what frequency range actually carries informa- 
tion on the reheating temperature. In other words, how wide of a band width is necessary 
to detect the knee shape with a good accuracy. We study this issue in Fig. 121 in which errors 
in Tr are plotted as a function of the upper frequency limit in the calculation of the Fisher 
matrix (/max in Eq- ©). Apparently frequencies above / ^ 0.3 Hz do not contribute to 
detection of the reheating temperature. This is because both the suppression of the signal 
amplitude due to reheating and the increase of the noise spectrum intensity prevents us from 
reaching the spectrum information. Thus, a moderate kink around is enough to reveal 
the reheating signature. 

Note that our analysis is performed imposing the consistency relation, r = — Sn^. Since 
the tilt of the spectrum is defined at the CMB scale, it can largely affect the amplitude at 
the direct detection scale as seen in Eq. ([9]). Therefore, an additional degeneracy between r 
and tlt arises and causes larger uncertainties in parameters if we take ut as a free parameter. 
For instance, in the case of r = 0.1 and Tr = lO^GeV with BBO/FP-DECIGO for 3-year 
observation, the uncertainty on the reheating temperature is degraded to = 2.9 x 10^ 
GeV by a factor of 8 compared to = 3.7 x 10^ GeV with fixed nx- However, combining 
with CMB B-mode polarization constraints would help to relax the degradation 



So far, we have fixed fiducial values of r and Tr. In the following, we discuss the fiducial- 
value dependence and predict the parameter space where the signature of reheating can be 
successfully detected. Fig. H] shows dependence on fiducial value of r. The marginalized 
error in Tr (cttr) is calculated by changing r with the fixed value of Tr = lO^GeV. The 
error becomes smaller as the gravitational wave background is detected with larger signal- 
to-noise ratio. According to Eq. ([9]), the amplitude of the spectrum at the frequency band 
of the experiment / ~ O.lHz takes the maximum value when r ~ 0.2 in the balance between 
the factors r (which increases the amplitude) and exp[— r ln(fc/A;o)/8] (which decreases the 
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FIG. 2: The 2a confidence level contours in the Tr — r plane for 1-year (dotted), 3-year (dashed) 
and 10- year (solid) observation by BBO/FP-DECIGO. The fiducial parameters are set as r = 0.1 
and Tr = lO'^GeV, which is shown by a cross mark. 

amplitude). This results in the smallest error on Tr around r ^ 0.2 in Fig. HJ Notice that the 
spikes originate from the fact that the Fisher matrix J>r oc {dQcwi/drY oc {1 — r ln(A;/A;o)/8} 
goes to zero around r ^ 0.23. This is an artificial effect due to our choice of parametrization. 

Similarly in Fig. |5]we show dependence on the fiducial value of Tr with the fixed value of 
r = 0.1. The error becomes smaller when the signature of reheating comes into the range of 
the sensitivity, which corresponds to the reheating temperature of about 10^ GeV to 10^ GeV. 
We should comment on our results at high Tr where the resultant ctt^ is underestimated. 
This is because the term of — 0.32a; ij in our approximated transfer function, Eq. f lT2|) . causes 
a small but an artificial "bump" in the spectrum, which does not arise in the full numerical 
calculation. For this reason, we highlight such questionable results at Tr > 10^ GeV as 
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r=0.1, TR=10^GeV (fR=0.26Hz) 




FIG. 3: The la marginalized errors are shown as a function of /max) calculated assuming r = 0.1 
and Tr = lO^GeV (corresponding to //? = 0.26 Hz) with 3-year observations. For clarity, the 
best-sensitivity frequency is plotted as a vertical dotted line. 



gray-shaded regions. The right panel of Fig. [5] shows a very promising fact in the case of 
Ultimate-DECIGO that the reheating temperature could be determined with 1% accuracy if 
1.2 X 10^ GeV < Tr < 3.3 x 10* GeV ^ for r = 0.1, and if 2.1 x lO^GeV < Tr < 7.0 x lO^GeV 
for r = 0.01. 

Finally in Fig. |6l we present the parameter space of r and Tr where the reheating 
signature is detected at greater than 2a level (i.e., T^/ctt^ > 2) with 3- year observation. 
Similarly to Fig. [5|, the region affected by our incorrect transfer function are indicated by 
the gray shade, and may shrink in a more realistic case. We also show the parameter region 
for the detection with a signal-to-noise ratio higher than 5 (S/N> 5) that is the criterion 



adopted in the previous work ISj. We cannot instantly compare two areas in the sense that 
the areas largely depend on variant criteria. However, the shape of parameter regions are at 
least different, originating from the fact parameter degeneracies are neglected in the previous 
work. We conclude that the parameter space presented here are more realistic and worth 



^ The upper value may be smaller by a several factor since the value is evaluated in the shaded region of 
the figure. 
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0.001 0.01 0.1 1 0.001 0.01 0.1 1 

r r 

FIG. 4: The marginalized la uncertainty in Tr as a function of r for BBO/FP-DECIGO (left 
panel) and Ultimate-DECIGO (right panel). The fiducial value of the reheating temperature is 
fixed to be Tr = lO^GeV. The upper horizontal axis represents the values of r^cw corresponded 
to r by Eq. dS]). 

while being pursued for detection of the reheating temperature with future gravitational 
wave detectors. 

IV. PROBING THE EQUATION OF STATE OF THE EARLY UNIVERSE 

In the previous section, we have focused on the determination of the reheating temper- 
ature by observing the change of the frequency dependence of the spectrum. On another 
front, direct detection of the inflationary gravitational wave background can be used to de- 
termine the equation of state of the very early Universe, at the cosmic temperature of around 
10^ GeV. So we could ask how sensitive the detectors are to different equation-of-state pa- 
rameters. Here, we calculate the Fisher matrix assuming that the frequency dependence of 
the spectrum is uniform over the range of sensitivity, and investigate how accurately direct 
detection can determine the equation of state of the early Universe. 
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FIG. 5: The marginalized la uncertainty in Tr as a function of Tr for BBO/FP-DECIGO (left 
panel) and Ultimate-DECIGO (right panel) . The fiducial value of the tensor-to-scalar ratio is fixed 
to be r = 0.1. 



DBCIGO Ultimate-DECIGO 



r 




Tr [GeV] Tr [GeV] 



FIG. 6: 2(7 detection region of Tr is shown as a blue shaded region for 3-year observations by 
BBO/FP-DECIGO (left panel) and Ultimate-DECIGO (right panel). The gray area represents 
the region where may be underestimated. In the light blue shaded region, the inflationary 
gravitational wave background would be detected with a signal to noise ratio higher than 5. 
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A. The spectrum with a generic equation of state 



Let us first relate the equation of state of the Universe w = p/p to the tilt of the 
gravitational wave background spectrum. If the initial power spectrum is assumed to have 
no tilt, A^ pj,;jj^ oc k^, the frequency dependence of the spectrum is determined only by the 
transfer function. Hence, Eq. ([8]) implies ^gw oc k'^T^{k). Since a gravitational wave, 
which has an initial amplitude of p^^^, maintains constant amplitude outside the horizon 
and it starts to decrease inversely proportional to the scale factor when the mode enters the 
horizon, the transfer function can be written as T^lk) = | /?.k,o I / 1 ^k,prim I = ('^o/'^in)~^; which 
means Qqw oc k'^af^. Combining the facts that the mode enter the horizon when k = aH 
and the Hubble expansion rate is given in terms of the equation of state as oc a^^'^^'^^\ 
we obtain oc /;;~2/(i+3«))_ 'jj-^^s, for modes which enter the horizon when the Universe has 
the equation of state w, the spectrum has the frequency dependence of 

2(3m-l) 

ficw oc k 3"-+! . (14) 

We parametrize the amplitude of the gravitational wave background spectrum, normalizing 

at / = F, as 

2{iw-l) 

^Gwif) = ncwAf/F)^^. (15) 



B. Result 

The fisher matrix is calculated by substituting Eq. (fT5|) into Eq. (|5]), taking i^cw.F 
and w as parameters. We investigate three fiducial models of the very early Universe; 
matter-dominated [w = 0), radiation-dominated {w = 1/3) and kination-dominated {w = 
1), which correspond to the frequency dependence of /~^, /° and respectively. Here, 
the normalization is taken at F = 0.203 Hz for BBO/FP-DECIGO and F = 0.158 Hz for 
Ultimate-DECIGO, which is chosen as the covariance matrix to be diagonalized for the flat 



spectrum 28| and is almost in the middle of the sensitivity range. Fig. [7] is one example of 
the error contours with the amplitude chosen to be f2Gw,F = 1-84 x 10~^^ corresponding to 
r = 0.1. From the figure, we see the constraint on w is significantly better in the matter- 
dominated case. This is because the measurable quantity in the experiment is not w but 
the tilt of the spectrum, tit = 2{3w — 1)/ (3w + 1). In this parametrization, the value of the 
tilt is more sensitive to the change of w when the fiducial model is w ~ than when w has 
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FIG. 7: The 2a confidence level contours in the w—Qq^^^f plane for 1-year (dotted), 3-year (dashed) 
and 10-year (solid) observation by BBO/FP-DECIGO. The fiducial value of f2GW,F is taken to be 
1.84 X 10^^® which corresponds to r = 0.1. The three panels represent a different model for the 
equation of state; w = (matter dominant, left panel), w = 1/3 (radiation dominant, middle 
panel) and w = 1 (kination dominant, right panel). These three fiducial points are shown as cross 
marks in each panel. 

larger value. This makes it easier for direct detection experiments to distinguish the model 
with smaller value of w. 

Fig. M shows dependence on the fiducial value of ^gw.f- The marginalized error 
is calculated by changing ^gw.f assuming the matter-dominated, radiation-dominated and 
kination-dominated universe, respectively. Note that the vertical axis does not represent 
0"^, but represents w. The blue shaded regions represent the la error ranges for the de- 
termination of w = 0,1/3 and 1, calculated assuming 1, 3 and 10- year observation with 
BBO/DECIGO and Ultimate-DECIGO. Obviously, the errors become smaller as the am- 
plitude of the gravitational wave background increases, which enables us to distinguish the 
value of w from that of the other models. 

V. CONCLUSIONS AND DISCUSSION 

In this paper, we study prospects for direct determination of the reheating temperature 
after inflation and the equation of state of the early Universe with future gravitational wave 
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FIG. 8: Dependence of the errors in w on the fiducial value of the gravitational wave amplitude, 
f^GW,F- The upper horizontal axis represents the values of r corresponded to f^GW,F by Eq. ([8]). 
The solid straight lines represent values of w for the matter (w = 0), radiation {w = 1/3) and 
kination (w = l)-dominated universe. The blue shaded regions indicate the la error range on each 
value of w for 1,3 10-year observations. Each three panels show the case where the fiducial model 
assumes the matter (left), radiation (middle) and kination (right )-dominated universe. The upper 
panel is for BBO/FP-DECIGO and the lower panel is for Ultimate-DECIGO. 

detectors. The reheating temperature, Tr, will be determined accurately for Tr ~ lO^GeV 
by BBO/FP-DECIGO, if the tensor-to-scalar ratio r is larger than ~ 0.05. Therefore, if the 
CMB B-mode polarization is measured by the Planck satellite or ground-based telescopes 
at the level of r ~ 0.1, we will have a good chance to detect the inflationary gravita- 
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tional waves at future space laser interferometers and determine/constrain the reheating 
temperature. Since Tr of 10^~^GeV is close to the upper bound from the gravitino problem 
in supergravity IstI], the detection of the gravitational waves will have an impact on the 
supersymmetric models if the Large-Hadron Collider finds signals of supersymmetry. De- 
termination of the reheating temperature is also closely related to the origin of the ordinary 
matter in the present Universe, since the amount of baryon asymmetry crucially depends 
on Tr for many baryogenesis mechanisms. For example, thermal leptogenesis scenario needs 



Tr > lO^GeV 



38 



39 



40|. 



and nonthermal leptogenesis needs Tr > 10^ GeV 
Moreover, future gravitational wave detectors have a potential to determine equation of 
state of the early Universe to a good accuracy. The matter, radiation and kination dominated 
universe are clearly distinguished by BBO/FP-DECIGO for r ~ 0.1. It will be a important 
step toward the complete understanding of the whole history of the Universe and the nature 
of the inflaton itself. 

As a final remark, we mention another probe of reheating. Accurate measurements of 



the CMB anisotropy are also sensitive to the reheating temperature 41 



42| . The expansion 



history of the Universe can be connected to the length of infiation and slightly affects the 
values of the slow-roll parameters. This is reflected both in the primordial spectra of the 
scalar and tensor perturbations, and is also applicable to direct detection of the inflationary 
gravitational wave background. Hence, direct detection helps to constrain the reheating tem- 
perature even in the case where the reheating signature is not seen in the frequency coverage 
of experiments 43|. Also, if there is a process called preheating 44 1, gravitational waves 
generated during preheating can be another probe of reheating 45H47j. These gravitational 
waves have a peak at the characteristic frequency, corresponding to the comoving Hubble 
scale at the end of inflation. Thus, they may provide us with complementary information 
on properties of inflation and reheating. 
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